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ABSTRACT 

Gamma-ray bursts (GRBs) are the most violent explosions in the Universe 
and can be used to explore the properties of high-redshift universe. It is believed 
that the long GRBs are associated with the deaths of massive stars. So it is 
possible to use GRBs to investigate the star formation rate (SFR). In this pa¬ 
per, we use Lynden-Bell’s c~ method to study the luminosity function and rate of 
Swift long GRBs without any assumptions. We find that the luminosity of GRBs 
evolves with redshift as L(z) oc g(z) = (1 + z) k with k = 2.43lo;3g. After correct¬ 
ing the redshift evolution through L 0 (z) = L(z)/g(z), the luminosity function 
can be expressed as ip(L 0 ) oc Lq °' 14±0 ' 02 for dim GRBs and ip(L 0 ) oc Lq °' 70±0 ' 03 
for bright GRBs, with the break point Lq = 1.43 x 10 51 erg s -1 . We also find that 
the formation rate of GRBs is almost constant at z < 1.0 for the first time, which 
is remarkably different from the SFR. At z > 1.0, the formation rate of GRB 
is consistent with the SFR. Our results are dramatically different from previous 
studies. Some possible reasons for this low-redshift excess are discussed. We 
also test the robustness of our results with Monte Carlo simulations. The distri¬ 
butions of mock data (i.e., luminosity-redshift distribution, luminosity function, 
cumulative distribution and log N — log S distribution) are in good agreement 
with the observations. Besides, we also find that there are remarkable difference 
between the mock data and the observations if long GRB are unbiased tracers of 
SFR at z < 1.0. 

Subject headings: gamma-ray burst: general - stars: formation - stars: luminosity 
function 
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Introduction 


Gamma-ray bursts (GRBs) are a kind of the most violent explosions in the Uni verse 


which radia t e huge energy in ga mma-ray in a short time (for reviews, see iMeszarosI 12006 
Zhang 2007; Gehrels et ah 2009). They are so bright and can be detected at much higher 


redshifts than supernov ae (SNe). Hitherto, the farthest GRB with spectroscopic redshift is 


GRB 090423 at z ~ 8.2 fTanvir et al. 


Wang 20131) , reionization ijTotani et a 


2009 


Salvaterra et al. 


(iWang et al.ll2012t ICastro-Tirado et al 


2006 


20131 


tools to probe the early universe ilBrorn rn fc Loebll2012t Wang et al. II20141) . i ncluding dark 
energy (IDai et al.ll2004 ISchaefeill20 07: W ang et al. 11201 lh . star form ation rate (lTotanilll997l; 


2009 a|). So GR Bs may be useful 


McQninn et ah 2008). and metal enrichment 


Theoretically, the progenitors of long GRBs with duration time Tgo > 2 s are thought 
to be collapse of massive stars (IWooslevI 19931). Observations also show that some long 


GRBs are associated with the deaths o 


massive stars which will give core collapse super- 


novae (Stanek et al. 

2003; 

Hiorth et al. 

20031). 

So GRBs can be used to 

investigate the star 

formation rate (SFR) at ’ 

ligh redshifts ( 

Totani 

1997; 

Wiiers et al. 

1998 

Lamb & Rcichart 

2000: Porciani & Madau 

2001; Bromrn & Loet 

2002 

; Lin e 

; al. 

2004; 

xistler et al. 

2009; 

Wang & Dai 2009; 

Wanderman & Piran 2010; 

Butler et al. 

2010; 

Elliott et al. 2012. 

2014). 


In order to measure SFR by using GRBs, the relation between the rate of GRBs and 
SFR should be known. Some studies found that the GRB rate is consistent with SFR 
at ab out z < 4.0, but has an excess at high redshift comparing with that expected from 
SFR fjLe &; Dermer 1 120071: iKistler et al. 1120091) . Some models have been proposed to explain 
the discrepancy betw een SFR and GRB rate, such as t he cosmic metallic ity evolution flLi 
20081: iQin et al.ll2010l) . superc onducting cosmic strings flCheng et al.l 120101) . evolving initial 


mass function (IMF 


break (Virgili et aL 


of s 


ars (Wang & Dai 20111). evolution of the GRB luminosity function 


2 01 lh. From the redshift distribution of GRBs and the metallicity of 


GRB host galaxies, IWang fe Dai (2 01 4) showed that the discrepancy between the GRB rate 
and SFR can be reconciled by considering that GRBs occur in low-metallicity galaxies. 


— ~ — ~ ^ ~ —--— -- 

nosity function and formation rate of GRBs (Fcnimore & Ramire Ruiz 

2000; 

Firmani et al. 

2004; 

Guetta et a 

1. 20 

05: Guetta & Piran 2007: Salvaterra & Chincarini 

2007; 

Salvaterra et al. 

2009b 

: Gao et al. 

2011 

). But the log N — log P distribution is the production of the intrinsic 


luminosity function convolving with the formati on rate of GRBs, so t he luminosity function 
and formation rate of GRBs are degenerated (IFirma/ni et al.l 120041: iGuetta fe Piranl 120071: 


Ho we ll et ah 2014). 


More ov er, there are several selection effects on the observed redshift distribution of 
GRBs (IC owardl 120071) . and so the rate of GRBs. The most important one is the observa- 
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tional limit of satellite. The Swift satellite has a flux limit, which means that it can not detect 
GRB dimmer than the flux limit. So the observed data is truncated, it will be difficult to get 
the intrinsic distribution of GRB before the selection effect is corrected. Lynden Bell (1971) 
applied a novel method to study the luminosity function and dens i ty ev olution from flux-limit 
quasar sample, which is called Lynden-Bcll’s c~ method flLvnden-Belllll971i). After that, this 
method has been w i dely used in some other objects, sue 


Merh 

illi et al. 

1986: 

Peterson et a 

2004|; 

Lloyd-Ronning et al. 

2002; 


1986 


i as galaxies 


Loti & Spillar 


Kocevski .V Liang 


(Kirshncr 


et ah 


1978 


1986). long GRB s (jYonetoku et al. 


2006 


Wu et al.ll20121) . and short GRBs 


(jYonetoku et al.l I2014J) . The basis of Lynden- Bell’s c~ method i s that the distributions o f 
luminosity L and redshift z are independent (jLvnden-Belll 11971c lEfron fe Petrosian Il992l) . 
So before applying this method, we need to test t he independence of L and z with a non- 
parametric test method provided by Efron & Petrosian 019921 ). Lynden-Bcll’s c~ method 
is a powerful method to estimate t he luminosity function and formation rate of objects 
with truncated sample. For example, Y onetoku et a l. (1200 4) used this method to derive the 
luminosity function and the formation rate of GRB from 689 BATSE GRBs with pseudo 
redshifts. They found that the GRB formation rate increases quickly at 0 < z < 1, and 
remains ap p roximately constant up to z ~ 10, which is consistent with observed SFR at 
z < 4.0. iKocevski fc Liang (120061) found that the GRB comovi ng rate density r ises steeply 
at z < 1.0, followed by flattening and declines at about z > 3.0. Wu et al. (2012) also stud¬ 
ied the formation rate of 95 Swift GRBs using Lynden-Bcll’s c~ method, and found that the 
GRB formation rate increases quickly at 0 < z < 1.0, and remains approximately constant 


, v ^ ~ __ 

Hopkins & Beacon 

2006; 

Kistlcr et al. 

2009; 

_ — 

Yiiksel et al. 

2008; 

Wang & Dai 

2009) 


In this paper, we study luminosity function and formation rate of latest Swift GRBs 
by using Lynden-Bell’s cr method. The Lynden-Bell’s c” method can break the degeneracy 
between luminosity function and GRB formation rate. This paper is organized as follows. 
We introduce the data from Swift satellite, and make the K-correction in the next section. 
The introductions to Lynden-Bcirs c~ method and the nonparametric r statistical method 
are given in section 3. In section 4, we derive the luminosity function and formation rate 
of GRBs. Monte Carlo simulation is used to test our results in section 5. Finally, section 
6 presents conclusions and discussions. Throughout the paper, we assume a flat ACDM 
universe with Q m = 0.27 and Hq = 70 km s _1 Mpc _1 . 
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2. GRB sample 


Swift is a multi-wavelength satellite to observe GRBs. It has instruments designed 
to analyze the bursts, X-ray and UV/optical afterglows, which can locate the positions of 
GRBs. We collect 127 long GRBs with well measured spectral parameters given by Fermi- 
GBM and Konus-wind. These GRBs have redshift data observed by Swift. In Table 1, we 
list the GRB sample, including name (Column 1), redshift (Column 2), low-energy power- 
law index a (Column 3), high-energy power-law index /3 (Column 4), peak energy of the vF u 
spectrum in observer’s frame (Column 5), peak flux in a certain energy range (Column 6), 
energy range (Column 7), bolometric luminosity (Column 8) and references (Column 9) of 
GRBs. 

Two spectral models are used to fit the spectra of GRBs, includi ng a p ower law with an 
exponential cutoff model (PLEXP) and the Band model (Band et al. 1993). The functional 
forms are as follows, 


m = a( 


100 keV 


(2 + a)E 

) a exp £ p , 


( 1 ) 


m = 



100 keV > 
E 


exp 
1 exp- 


(2 +a)E 


3—af { a —0)Ep 'jOt—ff 




( 2 ) 


100 keW v (2+a)100 keV' 

which represent a power law with an exponential cutoff model and the Band model respec¬ 
tively. 

Because the peak fluxes are observed in a lager range of redshifts which correspond to 
different range of energy bands in the rest frame_of JdRBs. The K-correction is required to 
get the bolometric luminosity of GRBs ( Bloom et al. 200ll ). The bolometric luminosity of 
GRB is 

■ (z)FK, (3) 


L = 47t d 2 L { 


where 


d L {z) = 


dz 


( 4 ) 


H 0 Jo \/l — + f^m(l + Z ) 3 

is the luminosity distance at redshift z, F is the peak flux observed between between a 
certain energy range (E min , E max ), and the K is the factor of K-correction. If the flux F is 
in units of erg cm -2 s -1 , the parameter K is defined as 


f E J nax Ef(E)dE 

J Hmin v 7 


K = 


( 5 ) 
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If the flux F is in units of photons cm 2 s x , the parameter K is defined as 


K = 


jYv/I [ { F; ] Ef(E)dE 


jf::HE)dE 


( 6 ) 


where f(E) is spectral model of GRB. Then we can obtain the bolometric luminosity L 
of each GRB. In Figure [U the blue dots show the bolometric luminosity of GRBs, and 
the line is the observational limit of Swift. The limit is cho sen a s a minimum flux F min = 
2.0 x 10 -8 erg s -1 cm -2 , which is consistent with that of Li (2008J). So the limit luminosity 
at redshift z is given as Li imit = 47 t d1(z)F m \ n . 


3. Lynden-Bell’s c method and nonparametric test method 


Lynden-Bcll’s c method is an efficient method to determinate the distribution of lumi¬ 
nosity and redshift of objects with truncated data sample, including quas ars (ILvnden-Bell 


1971 


1978 


2002 


Efron fe Petrosian 


1992 


Merighi et ah 1986 


Yonetoku et al 


Petrosian 11993 


Peterson et al. 1986 


2004 Wu et al 


2012 


Malonev &: Petrosian Ill999 ), g alaxies (Kirshner et al. 


Loh & Soilla: 


Yonetoku et al. 


19861) and GRBs (ILlovd-Ronning et al. 
2014 ). Lvnden-Belj f 1971 ) used 


this method to derive the luminosity function and density evolution from quasars with ob¬ 
servational selection for the first time. This method can break the degeneracy between 
luminosity function and formation rate. It is better to extract the luminosity evolution from 
the form of the luminosity function. If the parameters L and z in the distribution of lumi¬ 
nosity and redshift 4/(L, z ) are independent, we can rewrite ^(L,z) as ^(L,z) = ^(Vjf^z) 
(Efron & Petrosian 19921 ). where ip(L), the fraction of GRB brighter than L, is the cu¬ 
mulative luminosity function, and f>(z) is the redshift cumulative distribution. Unfortu- 
nately. the luminos i ty and the redshift of GRBs are not independent (lYonetoku et al. 1120041: 


Kocevski fe Liang 120061: IWu et al.l 120121) . so we should write the \I/(L, z) as Sif(L,z) = 


if z (L)f>(z) instead of 'f/(L,z) = if(L)f>(z ), where if z {L) is the luminosity function of GRB 
at redshift z. If we remove the effect of the luminosity evolution g(z), i.e. a transformation 
L 0 = L/g(z), the transformed luminosity L 0 are independent of redshift. As a result, we 
can obtain ^(L 0 ,z) = 4>(z)if(L 0 ). Using the relation L = L 0 g(z), we can write the ty(L,z) 
as ^(L,z) = ip z (L)(j)(z) = ijj(L 0 )^)(z). Then the luminosity function of GRB at redshift z is 
^z(L) = ip(L/g(z)). 


To get the form of luminosity evo lution g(z), we introduce a nonparametric test method 


proposed Ry 

Efron & Petrosian (1992). which has been widely used in previous literature 

(Petrosian 

1993 

; Ma 

onev & Petrosian 1999; Llovd-Ronning et al. 2002; Yonetoku et al. 

2004; 

Wu et al. 

2012|; 

Yonetoku et al. 

2014). For the ith data in the (L. z) data set, we can 














































































































define J, as fjEfron & Petrosian 1992) 

Ji = {j\Lj > Li, Z j < z™ ax }, (7) 

where L* is the ith GBR luminosity and zf iax is the maximum redshift at which a GRB with 
luminosity Li can be observed. This region is shown in Figure Q] as black line rectangle. The 
number of GRBs contained in this region is ri^ The number TV* = rq — 1, which means taking 
the '/’tli GRB out, is the same as c 


m 


Lvnden-Belll (1971). Similarly, J- is defined as 

,7'={j|L j >L“’V,< 2i }, (8) 

where L* im is the minimum observable luminosity at that redshift z*. This region is shown 
as red rectangle in Figure [H The number of events contained within this region is M*. 

We first consider the n, : GRBs in black rectangle in Figure [0 The number of events 
that have redshift z less than or equal to z; is defin ed as R ,. If L and z are independent, R* 
i s uniformly distributed b etween 1 and n t fjEfron fc Petrosian 1119921) . The test statistic r is 
( Efron fc Petrosian 1992 1 

(Ri - Ei) 


T = 




W, 


( 9 ) 


where E, = 


1 +rii 


Vi = 


nf-1 


are the expected mean and the variance of R, respectively. If 


2 ’ * 12 

the Rj is exactly uniformly distributed between 1 to n,, the samples of Ri < Ei and Ri > Ei 
should be nearly equal, and the test statistic r will nearly be 0. If we choose a from of g(z) 
that makes test statistic r = 0, the effect of the luminosity evolution can be removed by 
transformation of L 0 = L/g(z). 



dence between L 0 = L/g(z) and z by changing the value of k until the test statistic r is zero. 
Figure [2] shows the value of test statistic r as a function of k. From this figure, we find the 
best fit is k — 2.43l°;3g 1 a confidence level. So we take the luminosity evolution form g(z) 
as g(z) = (1 + z) 2 ' 43 . A hypothesis of no luminosity evolution k = 0, is rejected at about 
4.7cr confidence level. This value is similar with Yonetoku et ah ( 20041 ). whose k -value is 


k = 2.60ln'la hhd k = 0 is rejected with about 8.0cr significance from pseudo-redshift GRBs. 


Wu et ah (2012) also found that the value of k is k = 2.30lg; 


+0.56 
51- 


After converting the observed luminosity to non-evolving luminosity L 0 = L/(l + z) 2 ' 43 , 
we can derive the local cum ulative luminosity function L>(Ln) with nonp arametric method 
from the following equation f Lvnden- Belli 1 19 7li lEfron fe Petrosian" 1992 1. 
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where j < i means that GRB has luminosity Loj larger than L o*. The cumulative number 
distribution 0(z) can be obtained from 

m= ir+ jf)- (ii) 

j<i 3 


where j < i means that GRB has redshift Zj less than z t . The comoving differential form 
of 4>(z), which represents the cosmic formation rate of GRBs p(z), is more interesting. The 
formation rate of GRBs can be written as, 


P(z) 


d(j)(z) 

dz 


(! + *)( 


dV(z) 

dz 


-l 

1 


( 12 ) 


where (1 + z) results from the cosmological time dilation and dV(z)/dz is the differential 
comoving volume, which can be expressed as 


dv ( z ) _ 4vr (^)3( r dz ^ 1 

dz Hq Jo yR — f2 m 4 - I2 m (l + z)^ y/l — + fi m (l + ^) 3 


(13) 


4. Luminosity function and formation rate of GRBs 

In this section, we present results on the luminosity function and cosmic formation rate 
of GRBs. 


4.1. Luminosity function 


As discussed above, we get the form of luminosity evolution as g(z) = (1 + z) 2A3 by using 
the nonparametric r test method. The non-evolving luminosity Lq is defined as Lq = L/g(z), 
which is shown in Figure [3] Using this new data set, the luminosity function ip{Lo) can be 
derived by using Lynden Bell’s c~ method, which is shown in Figure 01 As shown in Figure 
01 the luminosity function ^j(Lq) can be fitted with a broken power law after removing the 
redshift evolution. The form of luminosity function i^(Lq) for dimmer and brighter bursts 
are fitted by 


V^o) oc 


T —0.14±0.02 

-^o 

r —0.70±0.03 
-^0 


U < Tq, 
Lq > Lq, 


(14) 


where L h n = 1.43 x 10 51 erg s 1 is the break point. This result is consistent with previous wor k 


fLlovd-Ronning et al. 

CM 

O 

o 

CM 

Yonetoku et al. 

2004; 

Kocevski & Liang 

2006; 

Wu et al. 

2012) 


It is necessary to point out that this luminosity function is only the present distribution 
























at z = 0 since the luminosity evolution is removed. The luminosity function ip z (L) at 
redshift z will be i^ z (L) = i\){L/g{z )) = + z) 2A3 ). So the break luminosity at 0 is 

L b z = L b 0 ( 1 + z) 2 - 43 . 


4.2. Formation rate of GRBs 


Figure 5 presents the cumulative GRB formation rate 4>(z). According to equation 
(fl2]l . in order to get the cosmic formation rate of GRBs, we need the differential form of 

cumulative number distribution dcj)(z)/dz. Figure El shows (1 + z)d(f)(z)/dz as function of 

redshift z. From this figure, we find that (1 + z)d(p(z)/dz increases quickly at z < 1, then 
keeps approximately constant for 1 < z < 4, and decreases sharply at z > 4 with a power-law 
form. But we are more interested in the comoving density rate. From equation f[T2jl . we can 
calculate the GRBs formation rate p(z), which is shown in Figure [71 In Figure [7] the blue 
stepwise line represents the comoving cosmic formation of GRBs as a function of redshift, 
and the error bar gives the la confidence level. The best-fitting power laws for different 
segments are 

f (1 + £)°- 04±0 - 94 z <\, 

p(z ) OC < (1 + ~j-0-94i0.il 1 < z < 4, (15) 

{ (1 + £)- 4 -36±0.48 z > 4 7 


with 95% confidence level. From this equation, we can derive the formation rate of GRBs 
at the local universe p(0) = 7.3 ±2.7 Gpc^ 3 yr -1 , which is larger th an prev i ous s tudies, 


e.g., p( 0) ~ 1.5 Gpc 3 
and p(0) > 0.5 Gpc -3 yr 


-i 


yr 

-l 


(Schmidt 

1999) 

p( o) 

(Pelaneeon et al 

2008) 


0.5 Gpc" 3 yr" 1 


(IGuetta et ah 20051), 


k eeps co nst an 


at low redshift in this paper, while it increases fast in other studies. But 

yr" 


Liang et al. (2007) found that the rate of low-luminosity GRBs is p(0) 


325 Gpc" 3 —- 1 


This local rate is not corrected for the jet beaming effect. 


Obviously, the formation rate o f GRB s is in contrast to previous e stima tes of the c omov- 
ing rate density by lYonetoku et al.l (120041) , iKocevski &; Liang (120061) and IWu et al.l (120121) 
with the same method. Our results shows that the formation rate of GRBs is almost con¬ 


stant at z < 1.0. But previous results give that 


z < 1.0 (lYonetoku et ah 2004; Kocevski fe L iang 


he formation rate in creases quickly at 


2006; IWu et al.l 120121 ) . which is consis¬ 


tent wi th SFR observat i on dHopkins fe Beacon 1 120061 1 . But our result is consistent with 


that of Petr osi a n et al. (1200911 well. Interestingly, the evolution of (1 + z)dcp(z)/dz shown 
in Figure El is consistent with the b ehavior of p(z) i n IWu et al. (1201211 . We als o test our 


Wu et al. (2012), and find 


program with the same GRB data of lYonetoku et al.l (120041 1 and _ 

that our (1 ± z)d(j)(z)/dz shows the same behaviors as the p(z) in lYonetoku et al. ( 2004 1 
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and IWu et al.l (120121) . So th ey might omit the d V(z)/dz term in their calculations. Be¬ 
sides, Yonetoku et ah (12004 ) and Wu et ah (2012) showed that the GRB rate increases as 
(T + z)® ' 0±1A an d (1 + z) 8,24=b4,48 at z < 1 respectively, which are much quickly than (1 + z) 2A 


of lKocevski fc Liang (120061) . 


Several comprehensive works studying the luminosity function and th e rate o f long 


GRBs have been recent ly published using different methods (such as IWanderman fc Piran 


2010 


Butler et al. 2QloK. A ssuming that the luminosity function is redshift independent, 


Wanderman fc Piran 


(120101 found that the power law index of luminosity function is 0.22 
at low luminosity, and 1.4 at high luminosity with break 10 52 ' 5 erg s^ 1 using long GRBs 
with redshifts determined from afterglow. The formation rate increas es as (1 + z) 21 up 


to z ~ 3.0 and it decreases as (1 + z)~ 1A at z > 3.0. 
the luminosity function is nearly flat oc L~ 

L~ 3 0 using a l arge sa mple of GRBs detected by Swift. The GRB rate is similar as that of 


Butler et ah (2010) derived that 
~ a2 below break 10 52 ' 7 erg s” 1 , and declines cx 


Wanderman fe Pi ran (120101 ). These results are different from our results. One reason is 


the GRB sample. We use the latest GRB sample, which have redshift observed by Swift 
and spectral parameters given by Fermi- GBM and Ko nus-win d. The luminosity function 


evolution may be the most important reason. IWanderman fe Piran 


Butler et al 


(20101) assumed no 


(hoiol ). no luminosity 


redshift evolution of luminosity function. In the fitting of 
evolution is required to produce the observed number of GRBs. But strong evolution of 
luminosity function is found in literatures . The evolution of lum inosit y can be parame terized 
as (l + z) L4 (Lloyd-Ronning et al. 2002). Yonetoku et al. ( 2004 ) and Wu et ah ( 2012) found 


that the evolution factor is g(z) = (1 + z) 2m and g(z) = (1 + z) 2,30 respectively. 


Tan et ah 


(120131) found that the luminosity function of GRB evolves with a redshift-dependent break 


luminosity L b = 1.2 x 10 51 (1 + z) 2 erg s 1 , 


Virgili et ah 


which is similar with our result. 

(20111) found that a evolution factor (1 + ^) 10±0 - 2 of luminosity function can fit the BATSE 
and Swift data. These works suggest that take a evolution factor into consideration is 
necessary. Besides, our GRB sample including those GRBs dimmer than 10 51 erg s’ 1 is 


another important reason. For example, Kistfer et ah ( 20081 ) found that the density of 
GRB is much higher at z < 1 if they included GRBs dimmer than 10 51 erg s 1 . In this 
work, we use Lynden-Bcll c~ method to correct the data truncated effect and consider the 
evolution of GRB luminosity function. So we don’t need to omit these dim GRBs. 

The relation between SFR and formation rate o f GRBs is attractive. We also compare 
our result with the observed SFR from|Hopkins & Beac on (120061 ) in Figure [8] Obviously, it is 
consistent with the observed SFR at z > 1.0, but remarkably different at z < 1.0. This trend 
means that the formation rate of GRBs p{z) does not trace SFR at low redshift z < 1.0. 
But at high redshif t z > 4.0, our result is consistent with the SFR d erived from GRBs 
( Yiiksel et al. 2008 : Wang &; Dai 2009 : Kistler et al. 2009 : Wang 2013 ). This result is dif- 
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ferent with others in previous literatures (jLlovd-Romming et ah 1120021: lYonetoku et al.l 12 004; 


Kocevski & Liang 

2006 

(2009 

). There are 

also 


Wu et al.l 120121) . but it is consistent with that of lPetrosian et al. 


tracers of SFR at low redshift. A stron g depe ndence of the GRB rate on host-galaxy prop¬ 
erties out to 2 ~ 1.0 is found bv IPerlev et al.l (120131). So use GRBs as dire ct tracers of the 
cosmic SFR is cautious at z < 1.0 (iPerlev et al.l 1201311 . Vergani et ah (2014) found that the 
mass distribution of long GRB host galaxy is different with the expected from star-forming 
galaxies observed in deep survey, which suggests that long GRBs are not unbiased tracers of 
star formation activity at least at z < 1.0. They also found that long GRB rate can directly 
trace the SFR starting from z ~ 4 and above. 


5. Testing with Monte Carlo simulation 

In this section, we use Monte Carlo simulation to test onr results. Firstly, we simulate a 
set of data (L 0 , z) which follows the distribution described by equation (1141) and equation (1141) 
using Monte Carlo method. Then, we transfer the luminosity Lq to L through L = Lo(l+z) k , 
where k = 2.43. So we can get sets of pseudo data of GRB luminosity and redshift (L, z). In 
the simulations, we create 200 pseudo samples. Each sample contains 130 GRBs. Then we 
use Lynden-Bell c~ method and nonparametric r test method to calculate the distributions 
of these pseudo samples. Finally, we compare the simulated data with observed data. 

Figure [9] shows the comparing results. The four panels give the luminosity-redshift 
distribution, luminosity function, cumulative distribution and log IV — log S' distribution. 
In the panel a, we randomly choose one pseudo sample of GRB from the 200 samples to 
compare with the observed data. The red dots and the bine dots represent the observed 
data and the simulated data, respectively. From this panel, we can see that the simulated 
data and the observed data have similar distributions. The other three panels b, c and d 
show the comparisons of the luminosity function, cumulative distribution and log N — log S 
distribution between the observed data and mock data. The red curves show the distributions 
of the observed data, blue curves give the distributions of all of the 200 pseudo samples of 
GRB data and the blue curves are the mean distributions of the 200 pseudo samples. We 
perform the Kolmogorov-Smirnov test between observed data and the mean distributions 
of simulated data. The chance probabilities of the three tests are 0.49, 0.86 and 0.96, 
respectively. From these panels, we can also conclude that the distribution of the observed 
data lie in the region of those pseudo data, which means that the derived luminosity function 
and formation rate of GRBs are correct. 

In order to test whether long GRBs are unbiased tracers of SFR at low redshift, we 
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simulate 200 n ew p seudo samples of GRBs by assuming that the GRB rate follows the SFR 


from Yuksel et al. (f2008j), i.e., p(z) oc (1 + z) 3A at z < 1, p(z) oc (1 + z) 0 3 at 1 < z < 4 


and p(z) oc (1 + £)~ 3 " 5 at z > 4. Then we use the same method to calculate the distributions 
of these pseudo data. We find that the cumulative redshift distribution of observed data 
is not consistent with the pseudo data, which is shown in Figure [TOl The red, blue and 
green curves have the same meanings as those in Figure [9j From Figure [TQ1 we can see that 
part of the cumulative redshift distribution line of observed data lies outside of the region 
occupied by pseudo GRB data, especially at z < 1.0. The Kolmogorov-Smirnov test between 
the distribution from observed data and the mean distribution of simulated data gives the 
chance probability of p — 6.9 x 10~ 12 . It means that long GRBs are not direct tracers of 
SFR at z < 1.0. 


Conclusions and Discussions 


In this paper, we use Lynden-BcH’s c~ method to study the luminosity function and 
formation rate of Swift long GRBs without any assumptions. First, we use a r statistical 
method to separate the luminosity evolution from the stable form of the luminosity function 


by choosing the evolution form g(z) = (1 + z) k . The most proper k is k = 2.43ln' 4 a, which 


gives r = 0. This v a lue is similar with those of lYonetoku et al.l (120041) . IWu et al.l (120121) and 
Kocevski &; Liang (12006 ). After correcting the luminosity evolution by L 0 = L/( 1 + z) 2A3 , 
the cumulative luminosity function ip(Lo) and cumulative number distribution <f>{z) of GRBs 
can be calculated, which are shown in Figure Q] and Figure [51 The luminosity function of 
GRBs can be well fitted with a broken power law form as ip(Lo) oc Lg°' 14±0 ' 02 and ip(Lo) oc 
for L 0 < L h 0 and L 0 > Lg respectively, where Lg = 1.43 x 10 51 erg s^ 1 is the break 


r —0.70=1:0.03 

point. 


We also derive the formation rate of GRBs through the differential form of the cu¬ 
mulative number distribution <j>{z). Figure [6] shows the evolution of (1 + z We find 
that (1 + z)^ft increases quickly at z < 1, then remains roughly constant at 1 < z < 4 
and finally decreases rapidly at high redshift. From equation (fT2l) . the cosmic formation 
rate of GRBs p(z) is derived, which is shown in Figure [3 The best-fitting power laws for 
different redshift segments are p(z) oc (1 + z) om , p(z) oc (1 + z)^ 0 ' 94 , p(z) oc (1 + z)~ 436 
for z < 1.0, 1.0 < z < 4.0 and z > 4.0 respectively. Our results show that the formation 
rate of GRBs is almost constant at z < 1.0. But prey i ous results give that the formation 


rate increa ses quickly at z < 1.0 (lYonetoku et al.ll2004 iKocevski fe Liang 120061: IWu et al. 


20121) . But lYonetoku et al.l (120041) and IKocevski fe Liang (120061) used the pseudo reds ihfts 
of GRBs rather than the observed redshifts. Besides, we find the p(z) in Wu et ah (2012) and 
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Yonetoku et ah ( 2004 ) increase fast at z < 1.0, which has the similar behavior of (1 + z)^^- 
shown in Figure El So they might omit the dV(z)/dz term in their calculations. From Figure 
IH1 it is easily to find that GRB formation rate p(z) is consistent with observed SFR at z > 1.0 
but entirely different at 2 < 1.0. It means that the formation rate of GRBs only traces SFR 


at z > 1.0, which is d ifferent with previous work (lYonetoku et al.l 12004 iKocevski fc Liang 


2006 1; Wu et aL- 2012). We find the low-redshift excess of GRB rate for the first time. 


Surprisingly, we find that formation rate of GRBs is consistent with SFR at z > 1.0, 
but shows an excess at low redshift z < 1.0 for the first time, which is different with previous 
works. Our result shows that formation rate of GRBs is larger than SFR at z < 1.0. Below, 
we will discuss some possible reasons for this low-redshift excess. 

The first one is that the definition of long GRBs is not c lear. In classical method, the 
long GRBs are defined by T 90 > 2 s (IKouveliotou et ah] 19931) . There is no clear boundary 


line in this diagram to separate the long and short GRBs. Moreover, Tg 0 is an observed 
time scale, which represents different time for GRBs at different redshifts. Meanwhile, the 


at z = 0.089, show no association of supernovae (IGal-Yam et al. 


Gal-Yam et al. 

2006; 

Fvnbo et al. 

2006; 


Gehrels et ah 2006). So more physical criterions are required to classify GRBs. B ecause 


only a subclass of GRBs can trace the SFR. Some attempts have been performed (IZhane 


20061: IZhang et al.l 120071.120091; iBloom et al.l l2008t iBromberg et al.l 120131: iLii et al.l 120141) . It 


has been suggested that GRBs can be classified physically into Type I (compact star origin) 
and Type II (massive star origin) (Zhang 2006; Zha ng e t al . 2007). 


The second one is that some selection effects have not been included in analysis. For 
example, it is easier to measure the redshift of those GRBs which are in lower redshift and 
therefore create a bias toward low redshift GRBs. It means that we lose some high redshift 
GRBs, so the formation rate of GRBs at low redshift we calculated will larger than the SFR. 
This bias can be removed by using samples with high completeness in the GRB redshift 


measurements. There are three of such a s ample in previous literature (IGreiner et al.l 12011 


Salvaterra et al. 

2012 ; 

Hiorth et al. 

2012 ) 


It could be considered in the future works. 


The third one is that there may exist a subclass GRBs, i.e., low-lumin osity GRBs 


(Cobb et al. 

2006; 

Pian et al. 

2006; 


Soderberg et al.l l2006t iLiang et al. 1120071). The local 


rate of low-luminosity G RBs may be high, i.e., p(0) = 100 — 1000 yr 1 Gpc 3 (jSoderberg et al. 


20061: ILiang et al. 1120071 ). much higher than high-luminosity GRBs. Th e progenitors of low 


luminosity GRBs may be different with those of high-luminosity GRBs ( Mazzali et al.ll2006 
Soderberg et al. 20061 ). The contamination from low-luminosity GRBs could lead to the 


low-redshift excess. 
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Table 1. List of long GRBs used in this paper. It gives the name, redshift z, spectra 
parameters a & f3, rest frame peak energy E p , peak flux F, energy range, bolometric 
luminosity L in 1 — 10 4 keV and reference of the parameters of spectrum of each GRB. 


GRB 

z 

a 


Ep(keV) 

Flux(erg/cm 2 / s) 
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L (erg/s) 

Ref. 
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050904 

6.29 

-i-i5±g:il 


314+2 

(1.84 ±0.41) X 10- 7 

15 - 

5000 

9 25”*" 2 ' 06 
y *^°-2.06 

X 

10 52 

4 

050922C 

2.198 

-0-83±g;U 


130.8+®:® 

4.5++ x 10- 6 

20 - 

2000 

1 95”*~°' 30 
i,yj -0.30 

X 

10 53 

1 

051001 

2.4296 

1 -i O+0.66 

0.56 


44-38++ 

L95 !o.57 x 10-8 

15 - 

350 

-1 qq+0.59 
-L-48-0.40 

X 

10 51 

5 

051109A 

2.346 

, 9 c+0.44 

i.ZQ-o 59 


16112 4 

5.8+ | x 10“ 7 

20 - 

500 

3-40++ 

X 

10 52 

6 

051111 

1.55 

-0-98+g 


179.70++ 76 

3-41++ x 10- 7 

15 - 

350 

7-04++6 

X 

10 51 

5 

060115 

3.53 

— 1 O"*” 0 ' 5 
1,u —0.5 


62+J 

0.9 ±0.1“ 

15 - 

150 

1 04+ 0,12 
1,u ^-0.12 

X 

10 52 

7 

060124 

2.296 

_i 29 -1-0 ' 14 

i.zy_ 0 ii 

9 9C+0.27 

^.z,o_o 88 

247.76+I+ 1 

2.66+; 7 ? x 10-6 

20 - 

2000 

-1 07+0.38 
1,0 ' —0.35 

X 

10 53 

8 

060206 

4.048 

— 1 12”*" 0 ' 3 

1,1 -0.3 


81+22 

° A -22 

(2.02 ± 0.13) x 10“ 7 

15 - 

150 

5 29 +0 ' 34 

O.z,»_o 34 

X 

10 52 

1 

060210 

3.91 

1 -| 9+0.26 

A * 1Z —0.26 


117+ 23 

AA/ -23 

2.8 ±0.3“ 

15 - 

150 

5.64+6° 

X 

10 52 

1 

060306 

3.5 

i 9+0.5 
_i - Z -0.5 


70+ 18 
' u -18 

(4.71 ± 0.278) x 10“ 7 

15 - 

150 

8-85++ 

X 

10 52 

1 

060428B 

0.348 

-0.94±i;gg 


21-7+4 

0.6 ±0.1“ 

15 - 

150 

9 1 r+0.36 
Z-±<J —0.36 

X 

10 49 

9 

060614 

0.125 

-l.57±g;“ 


302+ 44 

(4.5 ±0.7) x 10-6 

20 - 

2000 

9 qq+0.37 
0.79 

X 

10 50 

10 

060707 

3.425 

—0 66” 1 ” 0 ' 63 
U ' DO _0.63 


6612 

1.1 ±0.2“ 

15 - 

150 

1 1 9+O.2O 
A,AZ -0.20 

X 

10 52 

11 

060708 

1.92 

—0 93+ 0,47 

u.»o_ 0 43 


87.451“;“ 

1.78+;« x 10 “ 7 

15 - 

350 

5-S2+;? 7 

X 

10 51 

5 

060814 

0.84 

-1 43+ 0 - 16 
± -^°-0.16 


257+f 

(2.13 ±0.35) x 10-6 

20 - 

1000 

9-46+;gj 

X 

10 51 

12 

060908 

1.8836 

—0 93”*~°' 25 

u.»o_ 0 25 


14812 

(2.81 ± 0.23) x 10“ 7 

15 - 

150 

1 34+ 0 - 11 

X 

10 52 

1 

060927 

5.47 

n £1 +0-36 
U ' 8i -0.36 


71+ 14 

1 J-_i4 

(2.47 ± 0.17) x 10“ 7 

15 - 

150 

1 16+ 008 
1,10 _ 0.08 

X 

10 53 

1 

061007 

1.261 

—0 53+ 0 - 09 
u.oo_0 08 

—2 6i+°- 25 
z - o± -0.49 

498+? 

1 - 9 5+;li x 10-5 

20 - 

10000 

1 7Q+0-28 

-*-• *' °-0.22 

X 

10 53 

13 

061021 

0.3463 

i 99+O.I2 
x.zz-o i4 


77712? 

3.72+;?® x 10-6 

20 - 

2000 

1 76+ 0 - 25 
A -' °-0.77 

X 

10 51 

14 

061121 

1.314 

1 09+O.O4 
A-OZl —0.05 


606+° 

1.28+1? x lO" 5 

20 - 

5000 

1 48”*” 0-19 
-*-•^0-0.22 

X 

10 53 

15 

061222A 

2.088 

— 1 OO”*" 005 

±.UU_o 05 

9 09+O.38 
Z ' OZ -0.38 

353+ 

(4.8 ± 1.3) x lO” 6 

20 - 

10000 

1 48 - *" 0-40 

1 -^°-0.40 

X 

10 53 

1 

070110 

2.352 

_1 -t c+0.45 
i,10 -0.41 


108.33++? 2 

(5.168 ±0.831) x lO” 6 

15 - 

350 

2 95+ 0 - 87 

^.»o_o 59 

X 

10 51 

5 

070129 

2.3384 

I Oq+0.68 

-L.OO — 0 59 


65 96 +179 ' 79 
oo.ao_ 63 48 

2-72+® 55 x lO"® 

15 - 

350 

1 79+O.5I 
' -0.35 

X 

10 51 

5 

070306 

1.497 

—1.67±g+ 


> 105 

(3.04 ±0.164) x lO” 7 

15 - 

150 

> 1.04 x 10 52 

1 

070508 

0.82 

-o.8i±g ; g; 


1881® 

8.3++ x lO” 6 

20 - 

1000 

2 96”*"°' 37 
0.40 

X 

10 52 

16 

070714B 

0.92 

-0.86+1 


H20+ 780 

iizu_3 8 o 

2.7 ±0.2“ 

15 - 

150 

-1 99+O.O9 
A *^—0.09 

X 

10 52 

17,18 

070810A 

2.17 

—0 73+ 0 - 83 
u * /o _0.69 


42 23”*” 6 ' 62 

9.92l+ g x lO" 8 

15 - 

350 

4 92 +2 04 
^• yz -1.43 

X 

10 51 

19 

071003 

1.605 

—0 76+ 0,06 
u - 1 °-0.07 


780+J 

1.22+1? x 10- 5 

20 - 

4000 

9 1 0 + 0.34 
^• AO -0.39 

X 

10 53 

20 

071010B 

0.947 

9C+0.74 
i.zo-o 49 

—2 63+ 0 - 29 
z -°°-0.49 

52+4 

8.92++ x lO" 7 

20 - 

1000 

6 47”*" 2-17 
-4.34 

X 

10 51 

21 

071020 

2.145 

—0 65+ 0,27 
U.DO_o 32 


3221®? 

6.041+? x lO" 6 

20 - 

2000 

9 9C+0.45 
Z - ZO —1.44 

X 

10 53 

22 

071117 

1.331 

1 cq+0.15 
±.OO_ 0 16 


278+? 6 

6.661+? x lO" 6 

20 - 

1000 

9.95+;?? 

X 

10 52 

23 

071227 

0.383 

-0.7 


1000 

(3.5 ± 1.1) x lO" 6 

20 - 

1300 

9 ^9+0-79 

^■04 — 0 79 

X 

10 51 

24 

080207 

2.0858 

-1-I7±g;i; 

—3 59 +0 ' 45 

o.o»_o 45 

107-8++ 

1.0 ±0.3“ 

15 - 

150 

4 22”*” 1 ' 27 

^• zz -1.27 

X 

10 51 

25 

080319B 

0.937 

-0-86+2 

67512 

(2.26 ±0.21) x 10” 5 

20 - 

7000 

i o^ 0 - 10 

l.UO_Q 10 

X 

10 53 

1 

080319C 

1.95 

— 1 Q1+ 013 
A * UA -0.13 

1 q7+0.15 

-0.63 

307+? 1 

3.35++ x lO" 6 

20 - 

4000 

9 46”*” 2 ' 23 

X 

10 52 

26 

080411 

1.03 

! tr-,+0.04 
A * Oi -0.05 


259+? 

(1.28 ±0.16) x lO" 5 

20 - 

2000 

9-3311+ 

X 

10 52 

27 

080413A 

2.433 

—1 2 +01 
A - z -o.i 


170+ 80 

A ' u -40 

5.6 ±0.2“ 

15 - 

150 

4 41 +0.I6 
^ A -0.16 

X 

10 52 

28,29 

080413B 

1.1 

1 90+0.25 
A ‘^°—0.25 


78+? 

1.4 ±0.2“ 

15 - 

150 

1 ^^+0-06 
-L-^^-o.oe 

X 

10 52 

1 

080603B 

2.69 

-1.201+6 

9 rq+0.31 
Z.OO_o 84 

200+? 1 

1-51+.38 X lO" 6 

20 - 

1000 

1 11+0-29 
A ' 1A —0.28 

X 

10 53 

30 

080605 

1.6398 

-0-871+3 

297+? 

(1.6 ±0.33) x 10" 5 

20 - 

2000 

0 qq+0.69 
O.OO-o 69 

X 

10 53 

31 









Table 1—Continued 


GRB 

z 

a 


Ep(keV) 

Flux(erg/cm 2 / s) 

Range (keV) 

L (erg/s) 

Ref. 


080607 

3.036 

n «-+0.07 

—0.76_ o 06 

9 C7+O.I8 
-0.26 

348+ 7 

(2.69 ±0.54) x 10" 5 

20 - 4000 

2.21++ X 10 54 

32 

080721 

2.602 

—0 96 +0,07 
u.yo_o 07 

—2 42+0- 29 
-0.29 

4971“ 

(2.11 ±0.35) x 10- 5 

20 - 7000 

1.11++ x 10 54 

1 

080804 

2.2 

-0.88+1 


315.111+ 

(7.3 ± 0.88) x 10“ 7 

8 - 35000 

2.86++ x 10 52 

33 

080810 

3.35 


-2.5 

580+H 

1.71+ x 10- 5 

15 - 1000 

2-39++ x 10 54 

34 

080913 

6.7 

-0.821+® 

-2.5 

1 91 +232 
1^1_39 

(1.4 ±0.058) x 10"® 

15 - 150 

1.24++ x 10 53 

35,36 

080916A 

0.689 

—0 qq+ 0,95 
u.yy_ 0> 05 


208111 

(4.87 ± 0.27) x 10~ 7 

8 - 35000 

1 08++ x 10 51 

1 

081007 

0.5295 

_1 4+0.4 
A '^—0.4 


40+1 

2.2 ± 0.2° 

25 - 900 

4-35++ x 10®° 

1 

081008 

1.9685 

_f) oc+0.20 

U.OO_o 20 


176 4+ 23 - 9 

110.4±_ 23 9 

(3.21 ±0.33) x 10“ 7 

8 - 35000 

9-48++ x 10 51 

33 

081028 

3.038 

0 Qfi+0- 34 

U-OO—0.34 

-2-25+1 

59.66+jl 

(7.04 ±0.65) x 10 -7 

8 - 35000 

4.91++ x 10 52 

33 

081118 

2.58 

n RR+O-00 

—U.Do_ 0 09 


98.99+;gl 

(6.73 ±0.23) x 10“ 7 

8 - 35000 

3-99++ x 10 52 

33 

081121 

2.512 

—0 21“*" 0,28 
U-Z± —0.28 

-l.86±g;gg 

206.9++ 

5.1611H x 10" 8 

8 - 35000 

1.38+:“ X 10 53 

33 

081203A 

2.1 

_i 29”*” 0,15 


497+11 

3.71+jl x 10- 7 

15 - 350 

2.63l+| x 10 52 

1 

081222 

2.77 

—0 qn - *"® 93 
u.yu_ 0 03 

—2.33+;l 

167lf 

(1.76 ±0.058) x 10"® 

8 - 35000 

L01++ x 10 53 

1 

090102 

1.547 

—0 97+ 0 - 01 
u ' y '-0.01 


461+1 

(2.93 ±0.091) x lO" 6 

8 - 35000 

4.79++ x 10 52 

1 

090424 

0.544 

-1.02±g;gi 

_3 2fi”*” 0 ' 18 
o.zu_ 0 18 

1621+ 

(9.12 ±0.14) x 10“ 6 

8 - 35000 

1-14++ x 10 52 

1 

090429B 

9.4 

-0.69+$ 

46.21++ 6 

1.03++ x lO” 7 

15 - 350 

i.59+:I! x 10 s3 

19 

090516 

4.109 

—l.03±g;ff 

—21+ 01 
Z,1 -0.2 

51.4++ 

5.3 ±0.2° 

8 - 1000 

8-70+;!! X 10 52 

37 

090519 

3.85 

-0.5S±g;H 


120.5++ 

(2.25 ±0.17) x 10“ 7 

8 - 35000 

3.46++ x 10 52 

33 

090529 

2.625 

—0 715+O O 3 

u -' o -0.03 


199.9++ 

(3.006 ± 0.063) x 10“® 

8 - 35000 

1.82+;® 4 x 10 53 

33 

090618 

0.54 

—0 91+00 3 

uyi -0.03 

1 

to 

rf*. 

to 

1 + 
0 0 
b b 
-1 -a 

313.2+1° 

(1.73 ±0.073) x lO -5 

8 - 35000 

1-871+| x 10 52 

33 

090715B 

3 

1 1+0.37 
0.37 


134+1 

(9 ±2.5) x 10" 7 

20 - 2000 

8.781+1 x 10 52 

1 

090809 

2.737 

_n 47+O.O5 

u - 4 '-0.05 

t-1- 

q q 
b b 
+ 1 
CD 

CN 

1 

193.4+H 

(7.231 ± 0.6) x 10"® 

8 - 35000 

3.40l+| x 10 53 

33 

090812 

2.452 

— 1 03+0- 07 
i.uo_ 0 07 


586+“ 

2.77 ±0.28“ 

100 - 1000 

1.02++ x 10 53 

1 

090926B 

1.24 

-0.l9±g;gg 


95.61+ 

(4.73 ±0.28) x 10“ 7 

8 - 35000 

4.461+® x 10 51 

1 

090927 

1.37 

-0.68+$ 

9 1 O+0.01 
Z,1Z —0.01 

59.671+1 

(9.379 ± 0.23) x 10"® 

8 - 35000 

9.30I+! x 10 52 

33 

091018 

0.971 

1 c-q+0.48 
i.oo_ 0 .48 


28+1 

(4.32 ± 0.95) x 10“ 7 

20 - 1000 

4.90++ x 10 51 

1 

091020 

1.71 

— 1 20”* -0 ' 06 
1.ZU—006 

—2 29” 1 ” 0 ' 18 
z,.z,y_o 18 

187+1 

(1.88 ±0.026) x 10" 6 

8 - 35000 

3-441+ 4 x 10 52 

1 

091024 

1.092 

-1 5 +0 - 4 
± -°-0.4 

280+H 

3.461+1 x lO” 7 

15 - 350 

4-08+;“ x 10 51 

38 

091029 

2.752 

-l.46±g;i; 

9 1 9+0-01 
^•■ L ^-o.oi 

61.4ll+ 5 

1.8 ±0.1“ 

15 - 150 

1.40+;“ x 10 52 

39 

091127 

0.49 

-0-68+“ 

59.67+H 

(9.379 ± 0.23) x 10"® 

8 - 35000 

7.71++ x 10 51 

33 

091208B 

1.063 

— 1 29” 1 ” 0 ' 04 
i, ^ y -0.04 


119+ 

(2.56 ± 0.097) x lO"® 

8 - 35000 

1-811+6 x 10 52 

1 

100425A 

1.755 

-0-53++ 


< 36.02 

4-741+1 x lO" 8 

15 - 350 

< 1.41 x 10 51 

40 

100615A 

1.398 

— 1 24 -1 " 0 ' 08 
0.06 

9 97+O.II 
-0.12 

85 73+ 7 - 82 
oO. / 0—9 33 

8.3 ±0.2“ 

8 - 1000 

L061+! x 10 52 

41 

100621A 

0.542 

— 1 70 - *" 013 
±>/u -0.13 

—2 45+0- 15 
z, ^ o -0.15 

95+I 

(1.7 ±0.13) x 10“ 6 

20 - 2000 

3.24++ x 10 51 

1 

100728A 

1.567 

—0 47+0+ 5 

-0.15 

9 c+0.2 
z,o -0.3 

390+1 

(4.2 ±0.7) x lO”® 

20 - 10000 

6.45+;“ x 10 52 

42 

100728B 

2.106 

-0.90+“ 


130+ 

(5.43 ± 0.35) x 10“ 7 

8 - 35000 

1.97++ x 10 52 

1 

100814A 

1.44 

-0.55+3 

9 77+O.I7 

1 -0.17 

147+ 12 

A4W -10 

(7.5 ± 2.5) x 10~ 7 

20 - 2000 

l.° 8 l+6 X 10 52 

43 

100816A 

0.8049 

—0 31+0 05 
u-oi_ 0 05 

136.71+1 

15.59 ±0.25“ 

10 - 1000 

7.381+ 2 x 10 51 

44,45 

100906A 

1.727 

-1 1+ 0 - 1 
A,1 -o.i 

9 9+O.2 
^ -0.3 

180+1 

(2.7 ±0.3) x lO"® 

20 - 2000 

6.901+ 7 x 10 52 

46 

101213A 

0.414 

1 1+0.07 

A * — 0.07 

9 Qc+0.29 
0.72 

309.7++ 

4.67 ±0.32“ 

10 - 1000 

6.32++ x 10®° 

47 

101219B 

0.55 

r\ oq+0.36 
u.oo-o 36 

9 1 O+0.12 
0.12 

70+ 

2.0 ±0.2“ 

10 - 1000 

3.81++ x 10 50 

48 

110205A 

2.22 

1 c 9 +0.14 
A * oz —0.14 


222+4 

(5.1 ±0.7) x 10" 7 

20 - 1200 

2-651+! x 10 52 

1 

110213A 

1.46 

— 1 44+0.05 
a * 44 -0.05 

_ 2 f)5”*" 0 ' 28 

z -°°-0.62 

98.41+ 

17.7 ±0.5“ 

10 - 1000 

2-231+! x 10 52 

49 

110422A 

1.77 

-0.53+11 

246+1 

(1.2 ±0.15) x 10“® 

20 - 2000 

2.9°!!;!! x 10 51 

50 

110503A 

1.613 

-0-98+;“ 

9 ^+0.3 

' —0.3 

219+1 

(10 ± 1) x 10“® 

20 - 5000 

1-891!+ x 10 s3 

1 

110715A 

0.82 

1 90 + O.O9 
0.08 

9 7+O.2 
' —0.5 

120+i 

(1.1 ±0.1) x lO -5 

20 - 10000 

4.31++ x 10 52 

51 
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Table 1—Continued 


GRB 

z 

a 

p b 

E p (keV) 

Flux(erg/cm 2 / s) 

Range (keV) 

L (erg/s) 

Ref. 

110731A 

2.83 

n opo.03 

U -°—0.03 

-2-981“;® 

3041“ 

20.9 ±0.5“ 

10 -1000 

3 06 +0 ' 07 

O.UD—o 07 

X 

10 s ® 

52 

110801A 

1.858 

1 7+0.12 
' —0.15 

-2.5 

140l^ 70 

8.9ll?;f x 10" 8 

15 - 350 

4 43+ 0 - 84 

X 

10 51 

53 

110808A 

1.348 

— 1 07^” 0,12 
±,u '-0.11 


42381®“™ 

(1.1 ±0.2) x 10" 4 

20 - 10000 

8-96+!! 

X 

10 53 

54 

110818A 

3.36 

1 qq+0.08 

x.oo-o 08 


256.3l®®j 

5.0 ±1.4“ 

10 - 1000 

7 9q+2.03 

' - ZO -2.03 

X 

10 52 

55 

111008A 

4.9898 

— 1 36” 1 ” 0,24 
l.OD-o 21 


149188 

(1.4 ±0.3) x 10-® 

20 - 2000 

4 95+ 106 

X 

10 53 

56 

111107A 

2.893 

1 OQ + O.^l 

l.dO_o 21 


1081“ 

2.6 ±0.3“ 

10 - 1000 

1 SI +9-21 

X 

10 52 

57 

111123A 

3.1516 

— 1 30 +0,26 
l.OO-o 24 

9 7+0.3 
' —0.3 

107.79l“5 Q 38 

7-89l8;J 9 x 10- 8 

15 - 350 

9-8°±i;“ 

X 

10 51 

40 

111228A 

0.714 


34I3 

27 ±1“ 

10 - 1000 

6 67+ 0 - 25 
-0.25 

X 

10 51 

58 

120119A 

1.728 

-0.98±g;“ 

9 qc+0.09 
^.ou_o 09 

189.21®® 

16.86 ± 0.39“ 

10 - 1000 

5 98+ 014 

O.UO—o 14 

X 

10 52 

59 

120326A 

1.798 

-0.98±g;ii 

9 cq+0.15 

z.oo_ 0 15 

46.451®;®' 

3.1 ±0.05“ 

10 - 1000 

5 91+ 010 
°- y± -o.io 

X 

10 51 

60 

120327A 

2.813 

— 1 14+ 0 - 26 
J - ±f *-0.28 

1 s+0-^ 
i ‘ 8 -0.2 

106.09l®°® 6 

3.88lgJ| x 10- 7 

15 - 350 

0 49+ 0 - 57 

-0.46 

X 

10 52 

40 

120712A 

4.1745 

—0 6” 1 ” 0-2 

U.O-o 2 

124188 

3.5 ±0.2“ 

10 - 1000 

qc+0.08 

x.oo-o 08 

X 

10 53 

61 

120714B 

0.3984 

-0.29lpg 6 


nn 0+25.92 
t>U.O_io 22 

l.SlJll x 10-8 

15 - 350 

1 ik+0.72 
1 - 1O -0.44 

X 

10 49 

40 

120724A 

1.48 

- 0 . 75 +lH 


< 31.9 

1.951111 x 10-8 

15 - 350 

< 4.19 x 

1 0 50 

40 

120802A 

3.796 



52.96lg 2 g4® 

n85lo 27 x 10-7 

15 - 350 

q ki +U.6b 
°- O± -0.51 

X 

10 52 

40 

120811C 

2.671 

— 1 19+ 032 
x,xy -0.30 


46.26+“ 

2.47ll;?l x 10- 7 

15 - 350 

9 9q+0.23 

z.zO-o 20 

X 

10 52 

40 

120907A 

0.97 

—0 73+ 0 - 25 

U ‘' °-0.25 


154.51“;® 

4.3 ±0.4“ 

10 - 1000 

9 cq+0.24 
Z.OO_ 0 24 

X 

10 51 

62 

120909A 

3.93 

_1 3+ 01 

x -°-o.i 

9 q+0.1 

z -O-0.1 

370111° 
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Notes: a For these GRBs, the peak flux is in units of photons cm -2 s —1 . 

b For those GRBs with value, the spectra are described well by Band model. But for the GRBs without [3 value, the spectra 
are describe d by power- l ow wi th an exponential cutoff mod el. _ 
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Fig. 1.— The luminosity distribution of 127 GRBs after K-correction. The blue dots repre¬ 
sent GRBs and the blue line represents observational limit of Swift. We take the flux limit 
as 2.0 x 10 -8 erg cm~ 2 s' 1 . Mi and N t are also shown. The error bars are la errors. 
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Fig. 2.— The value of test statistic r as a function of k. The red dotted line represents the 
best fit for r = 0, and the black dotted lines are lcr errors. The value of k is k = 2.43lg 3 g 
lcr confidence level. It also shows that r = 4.7 when k = 0, which means k = 0 is excluded 
at 4.7cr confidence level. 
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Fig. 3.— The non-evolving luminosity L$ = L/{ 1 + z) 2A3 of 127 GRBs above the truncation 
line. The blue line represents observational limit. The error bars are lcr errors. 
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Fig. 4.— The cumulative luminosity function ijj(L 0 ), which is normalized to unity at the 
lowest luminosity. The red line is the best fit with a broken power law model. The luminosity 
function can be expressed as 4>(L 0 ) oc Lq°' 14±0 ' 02 for dim GRBs and ^(L 0 ) £” 0 ' 0±0 03 f or 

bright GRBs, with the break point L b 0 = 1.43 x 10 51 erg s _1 . 
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Fig. 5.— The normalized cumulative redshift distribution of GRBs. 
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Fig. 6.— The evolution of (1 + as function of redshift z with 1 a errors, which is 

normalized to unity at the first point. 
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Fig. 7.— The comoving formation rate p(z) of GRBs, which is normalized to unity at the 
first point. The lcr error is also shown. 














31 



Fig. 8.— The comparison between GRB formation rate p(z) (blue) and the observed SFR. 
The SFR data are taken from Hopkin s & Beacom (2006), which are shown as red dots. The 
SFR data from lBouwen s e t ah (20111) (stars) and IWang I (120131 ) (open circles) are also used. 


All error bars are la errors. 
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Fig. 9.— The comparison of the simulated data (blue) and the observed data (red). These 
four panels show luminosity-redshift distribution, cumulative luminosity function, cumulative 
number distribution and logiV — log S' distribution respectively. For panel (a), we choose 
one sample from the 200 simulated samples randomly. For other panels (b), (c) and (d), the 
green curves represent the mean distribution of those 200 simulated samples. The chance 
probabilities of Kolmogorov-Smirnov tests between the distributions of observed data and 
the mean distributions of the simulated data are also presented. 
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Fig. 10.— The comparison between the cumulative redshift distributions of simulated data 
(blue) and observed data (red). The mean distributio n of the 200 simulated samples is 


given by green curve. In this case we use the SFR from Yii ksel et ah ( 2008! ). The chance 


probability of the Kolmogorov-Smirnov test between the observed data and the simulated 
data is p — 6.9 x 10~ 12 , from which the observed data and the simulated data from the same 
sample can be discarded. 




























